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The LHCb collaboration has recently performed a first measurement of the angular production 
asymmetry in the distribution of beauty quarks and anti-quarks at a hadron collider. We calculate 
the corresponding standard model prediction for this asymmetry at fixed-order in perturbation 
theory. Our results show good agreement with the data, which is provided differentially for three 
bins in the invariant mass of the bb system. We also present similar predictions for both beauty-quark 
and charm-quark final states within the LHCb acceptance for a collision energy of -yi = ISTeV. 

We finally point out that a measurement of the ratio of the bb and cc cross sections may be useful 
for experimentally validating charm-tagging efficiencies. 


INTRODUCTION 

The LHCb collaboration has recently performed a first 
measurement of the angular asymmetry in bb produc¬ 
tion at a hadron collider [1]. More specifically, LHCb 
has measured the forward-central asymmetry of 6-quark 
pairs with lfb“^ of Run-I data, collected at a 

centre-of-mass energy (^/s) of 7TeV in pp collisions.^ 

Instrumented in the forward region, the LHCb detec¬ 
tor operates in a kinematic regime which is well suited 
to measure heavy-quark production asymmetries. For¬ 
wardly produced heavy-quark pairs, particularly at high 
invariant mass, are sensitive to colliding partons at both 
moderate and large momentum fractions within the pro¬ 
ton. This kinematic sensitivity provides a unique oppor¬ 
tunity to perform asymmetry measurements as the dilu¬ 
tion from the otherwise overwhelming symmetric gluon- 
gluon-fusion (gg) production mechanism is reduced to a 
manageable level. Moreover, the ability of the LHCb de¬ 
tector to efficiently tag semi-leptonic B decays has made 
this measurement possible with the available data set. 

A notable feature of is that it receives a large cor¬ 
rection from purely electroweak (EW) effects when the 
invariant mass of the 6-quark pairs is close to the Z- 
boson resonance, as pointed out in [3]. Measurements 
of ApQ are therefore not only of general importance as a 
test of the standard model (SM), but also provide valu¬ 
able model-building input 11HI2], serving to restrict the 
set of new-physics scenarios which were suggested as an 
explanation of the anomalously large forward-backward 
asymmetry in top-quark pair production as observed at 
the Tevatron [TSHTB] . Although the tensions between the 
experimentally observed top-quark asymmetries and the 
corresponding SM predictions have been to a large extent 
resolved, owing to experimental improvements m and 


^ Subsequently, the forward-backward asymmetry of 6-quark pairs 
has also been measured in pp collisions at the Tevatron [313]. 


the inclusion of next-to-next-to-leading order (NNLO) 
QCD corrections m, measurements of App, remain in¬ 
teresting in view of the persistent discrepancy between 
the experimental data and the SM predictions for the 
Z ^ bb pseudo-observables (TH]. 

In this letter we report on the SM calculation of Ap^, 
comparing our results to the LHCb data. We also provide 
predictions for Ap^ and the corresponding charm-quark 
pair production asymmetry Ap^ at ^/s = 13TeV, rele¬ 
vant for data taking commencing in the next two years. 
Our SM calculations improve on work very recently pre¬ 
sented in [T^ . 

ANATOMY OF ASYMMETRY 

The forward-central asymmetry in heavy-quark pro¬ 
duction is defined in terms of the pp QQ cross sec¬ 
tion a in the following way 

aQQ ^ O' (Ay > 0) - ct (A^ < 0) 
cr {Ay > 0) -k (T {Ay < 0) ’ 

where Ay = |j/q| — |?/q| is the difference of the absolute 
rapidities of the heavy quark Q and anti-quark Q, evalu¬ 
ated on an event-by-event basis. At the LHC, the pp ini¬ 
tial state is symmetric with respect to the direction of the 
incoming proton. However, an asymmetry is present in 
the momentum fraction distributions of quark and anti¬ 
quarks within the proton as a consequence of the valence 
content, and the definition of Ap^ is chosen to reflect this 
— i.e. to gain sensitivity to the underlying parton distri¬ 
bution function (PDF) asymmetry in the incoming quark 
momenta \pz {g) \ — \pz (?) | ■ Consequently, the asymmetric 
contribution to the numerator of 0 arises from subpro¬ 
cesses of the form gg -T QQX, where X denotes a partic¬ 
ular final state such as g or 7. Note that by crossing sym¬ 
metry, the subprocesses g{g)X QQg{g) also contribute 
to the numerator. As demonstrated first in [201 1^ - 
the dominant contribution to the top-quark asymme¬ 
try arises at next-to-leading order (NLO) in QCD. The 
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contributions from mixed QCD-electroweak corrections, 
which are structurally similar to the NLO QCD effects, 
and pure EW contributions, which arise at leading or¬ 
der (LO) due to the presence of both vector and axial- 
vector couplings in the subprocess qq ^jZ —>■ QQ, 
have also been computed In the case of tt pro¬ 

duction at LHCb, the mixed QCD-EW contributions 
amount to 0(10%) of the total QCD corrections, while 
pure LO EW effects are negligibly small. In contrast, 
both the bb and cc asymmetries can receive large contri¬ 
butions from pure EW effects when the invariant mass 
mgg of the heavy-quark pair is in the vicinity of the Z 
pole. In this kinematic regime, higher-order corrections 
to the LO EW contribution can in principle also become 
important, and should be included if possible. 

The prediction for is cast in terms of a Taylor 

series expansion in powers of the strong (a^) and the 
electromagnetic (a) couplings in the following way 
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Here the terms cTa^°^ and correspond to the asym¬ 

metric NLO QCD and mixed NLO QCD-EW contribu¬ 
tions, respectively, while and (Tq*'^^ represent the 

pure EW asymmetric contributions and the correspond¬ 
ing leading QCD correction. In the denominator, our 
calculations include the LO symmetric QCD and pure 
EW contributions and as well as the associ¬ 

ated QCD corrections and 

Analytic formulas for the term CTo^*^^ can be found 
in m- Approximate results for the contribution 
are also provided in that article, but these results are 
not applicable in the resonant region mgg ~ mz, which 
is relevant for the beauty-quark and the charm-quark 
asymmetries. We have therefore computed the relevant 
formula for with the help of FeynArts m and 

FormCalc [28]. The numerical integration of these for¬ 
mulas is performed using the Vegas algorithm as imple¬ 
mented in the Cuba library [29] . and the complex scalar 
one-loop integrals are evaluated with the OneLOop pack¬ 
age |3D|. The contributions and have been 
calculated utilising the helicity amplitudes of m which 
have been extended to include the 0{as) corrections as¬ 
sociated to the final-state heavy-quark lines. All of the 
aforementioned terms have been calculated with physi¬ 
cal heavy-quark masses, with the exception of the term 
(Tq;s which is computed in the massless limit mg = 0. 
The final contribution ^ is obtained with the matrix 
elements of |32|, that are incorporated in POWHEG |33|. 
To combine all of our predictions, we perform a change 
of renormalisation scheme for the symmetric QCD com¬ 
putation which is performed in a fixed-flavour-number 
scheme (the four-flavour scheme for 6-quark pair pro¬ 


duction). This is done following the procedure outlined 
in [34], allowing the contributions from all subprocesses 
for both beauty and charm predictions to be consistently 
convoluted with five-flavour PDFs. We note that sev¬ 
eral cross checks of our predictions were performed with 
MCFM [ 33 . Further details on our computations, includ¬ 
ing analytic formulas for all asymmetric terms will be 
presented elsewhere [36] . 

Our numerical results are obtained for the following 
choice of input parameters [37] • = 173.25 GeV, mb = 

4.75 GeV, = 1.5 GeV, mw = 80.385 GeV, Tw = 
2.085 GeV, mz = 91.1876 GeV, Tz = 2.4952 GeV and 
Gf = 1.16638 • 10“®GeV“^. To describe the Z reso¬ 
nance, we adopt the complex-mass scheme (see e.g. [35] ') 
and determine the sine of the weak mixing angle and 
the electromagnetic coupling from = 1 — m^/m^ 
and a = y^/TTGpm'^s^, respectively. All contributions 
to are computed with NNPDF2.3 NLO PDFs [39] us¬ 
ing as{mz) = 0.119. Finally, we evaluate the ratio in (H 
exactly as written, rather than performing a further ex¬ 
pansion of higher-order terms in the denominator. These 
are our most advanced theoretical predictions, and we 
will refer to the results obtained in this way as “NLO” 
asymmetry predictions. In addition, to explore the con¬ 
vergence of the perturbative series, we also compute the 
leading contribution to the asymmetry with NNPDF2.3 
LO PDFs and as{mz) = 0.119. This corresponds to 
dropping the terms (t|^^^ and (Ta[l^ in the expression (|^, 
and results obtained in this way are therefore referred to 
as “LO” asymmetry predictions. For both LO and NLO 
predictions, a scale uncertainty is evaluated by indepen¬ 
dently varying the factorisation and renormalisation 
/iij scales by a factor of two around m^, imposing the 
constraint 1/2 < < 2. The total uncertainty of 

the predictions for the asymmetry is then found from the 
envelope of the different results. 


COMPARISON WITH ^/s = 7TeV DATA 

To compare our predictions to the available data, a 
fixed-order analysis^ is performed which includes the ap¬ 
propriate experimental cuts to mimic the LHGb selec¬ 
tion requirements. Jets are clustered using the anti- 
kt algorithm [41] with a distance parameter R = 0.7. 
Since in the data the reconstructed jets are corrected 
to the parton level, this procedure should allow for a 
fairly good comparison. The reconstructed jets are re- 


^ We have investigated resumming potentially large corrections 
arising from the soft and small-mass logarithms associated with 
energetic heavy-quark production as calculated in m, but it is 
not clear that results obtained in the soft-gluon emission limit 
are applicable in the case considered here. In consequence, we 
do not include such effects in our analysis. 
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FIG. 1: NLO and LO predictions of the beauty-quark 
forward-central asymmetry at ^/s = 7TeV within the LHCb 
acceptance. The statistical and systematic uncertainties of 
the measurement have been added in quadrature to obtain 
the shown experimental error bars. For further details con¬ 
sult the text. 


quired to be within the pseudo-rapidity range 2 < ry < 4, 
to have a minimum transverse energy of Et > 20 GeV 
and are also constrained to have an opening angle of 
A(/) > 2.6 in the transverse plane. The calculation of 
is then performed differentially in the invariant mass bins 
G [40, 75] GeV, [75,105] GeV and > 105 GeV to 
match the LHGb analysis. 

Our NLO and LO predictions for the asymmetry are 
compared to the existing LHGb data in Fig. [^ and a 
numerical summary of these predictions is provided in 
Table [T] It is evident that the SM prediction and the 
data are in full agreement within the given uncertainties, 
which are currently dominated by the statistical preci¬ 
sion of the measurement. Besides this, there are several 
important features of the theoretical prediction which we 
would like to mention briefly. 

Similarly to the top-quark asymmetry, the dominant 
contribution to Hpp. arises from the QGD contribu¬ 
tion for most values of the invariant mass of the 
&-quark pair, the important exception being the central 
mass bin with G [75,105] GeV. In this region, the 
double-resonant contribution from Z-Z interference be¬ 
comes dominant, accounting for the bulk of the total 
asymmetry, as shown in Table |lj The QGD correction 
to this contribution in the resonant bin is observed to 
slightly decrease the numerator, while the impact on the 
denominator is marginal. The contribution from mixed 
QGD-EW and ^-Z contributions are numerically negligi¬ 
ble in this bin (in the case of the QGD-EW corrections 
this confirms the recent estimate of (T^] ), which is a con¬ 
sequence of integrating terms that have a single Z prop¬ 
agator over the resonant region. The size of the mixed 


NLO 

[%] 

QGD 

QCD-EW 

EW 

rrifji G [40, 75] GeV 

0.59«:f6 

100.6% 

-4.9% 

4.3% 

myf, G [75,105] GeV 

Q 00+0.09 
Z.ZO_o.75 

33.5% 

-1.4% 

67.9% 

rrif^i > 105 GeV 

-L.uy_Q 72 

86.6% 

-7.1% 

20.5% 

LO 





G [40, 75] GeV 

Q oc+0.04 

U.OD_o.o3 

105.0% 

-5.1% 

0.2% 

mt,i G [75,105] GeV 

0 oq+0.45 

Z.OO_Q 37 

30.9% 

-1.2% 

70.3% 

rrifji > 105 GeV 

1 Q4+0.12 

96.8% 

-8.3% 

11.5% 


TABLE I: NLO and LO predictions of the bb forward-central 
asymmetry at y/s = 7 TeV within the LHCb acceptance. The 
relative contributions to Apc from QCD, mixed QCD-EW 
and EW corrections are also provided. 

QGD-EW corrections is further reduced across all bins by 
a partial cancellation between contributions arising from 
uu and dd initial states, which is a result of the u and 
d quarks having opposite electromagnetic/weak charges. 
We also note that the gg-initiated contribution is numer¬ 
ically sizeable, accounting for up to almost 10% of the 
entire term. 

The coloured bands around the central values of 
the NLO (yellow) and LO (green) asymmetry predic¬ 
tions shown in Fig. are due to scale variation alone. 
Away from resonance, the scale uncertainty of the LO 
asymmetry is artificially small as both numerator and 
denominator of Ap^ are dominated by the leading QGD 
contributions and respectively. In both cases, 
the dependence enters only via the scale dependence 
of as, and two powers of as cancel in the ratio. For 
rrifji G [75,105] GeV, the LO scale uncertainties are more 
pronounced as the dominant correction in that bin 
has no fj-u dependence, while the denominator is to first 
approximation given by 

Given the unreliable evaluation of theoretical uncer¬ 
tainties in the case of the LO results, and the ambigu¬ 
ity that these results strongly depend on whether NLO 
or LO PDFs are used for the simultaneous computation 
of the numerator and denominator, we believe that the 
NLO predictions are most trustworthy. Let us finally 
add, that although our NLO and LO predictions are 
consistent within uncertainties, this feature arises due to 
the presence of a hard LO gluon PDF at large-a; which 
partially compensates the effect of missing higher-order 
corrections in the symmetric cross section |26j . Given 
that the gluon PDF is essentially decoupled from many 
observables in a LO global PDF fit, this compensation 
can be regarded as accidental. Conversely, if the LO 
predictions are obtained using NLO PDFs, like for in¬ 
stance in m, the resulting LO asymmetry predictions 
tend to be systematically higher by about 50%. Further 
improvement beyond the NLO predictions would require 
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0 0 
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TABLE II: NLO and LO predictions for the bb forward-central TABLE IV: Ratios of charm and beauty predictions based on 
asymmetry at y'i = 13TeV. Predictions for the symmetric LHCb acceptance cuts and = 13TeV. 
cross section within the LHCb acceptance are also given. 


the inclusion of 0{ag) corrections to both the numerator 
and denominator in ([^, as done in m for tt production. 
Such a computation is beyond the scope of this letter. 


PREDICTIONS FOR = 13 TeV 

In addition to bb asymmetry measurements, LHCb 
has recently demonstrated the ability to efficiently tag 
charmed jets |42j . which in the future may allow for 
asymmetry measurements involving cc final states, as¬ 
suming that charge tagging can be established for charm 
quarks. We therefore also provide predictions for at 
i/s = 13 TeV, relevant for data taking in the initial years 
of LHC Run-H. 

Following the y/s = 7 TeV analysis strategy, we provide 
our fixed-order predictions in bins of the heavy quark- 
pair invariant mass bins rngQ, and also apply the same 
Et > 20GeV, 2 < 77 < 4, and A(j) > 2.6 cuts. The total 
asymmetry and symmetric cross section predictions for 
the beauty final state are provided in Table [IT| The rel¬ 
ative contribution to the asymmetry from QCD, mixed 
QCD-EW and EW corrections is qualitatively unchanged 
with respect to the ^/s = 7 TeV predictions (see Table|^. 
Importantly, at = 13 TeV, the symmetric cross sec¬ 
tion increases by a factor of three to five, while the total 
asymmetry is reduced by approximately a factor of two. 
This implies that, with 5fb~^ of integrated luminosity, 
LHCb should be able to improve the statistical precision 
of their bb measurement by at least a factor of two. This 
assumption conservatively assumes no improvement in 
the tagging efficiency. 



[GeV] 

[40,75] 

[75,105] 

> 105 

NLO 

Af-1: [%] 

af [nb] 

0.40l°;?^ 

57.9lf5;° 

n '7Q+0.09 
^°-0.20 

13.6115 

0.89l5'i? 

4.0315:°® 

LO 

A^^ [%] 

0.201^2 

0.90151? 

n 70 + 0 .O 6 
'“'• ‘°-0.06 

aj" [nb] 

llb+m:3 

12.515:? 

4 1 c+0.98 

+.-LO_o.74 


TABLE III: Predictions for charm-quark pair final states ana¬ 
logue to Table [III 


Our NLO and LO predictions for cc production are 
summarised in Table HI As pointed out recently in 
there is an additional contribution arising at 0{asa) due 
to s-channel gluon exchange interfering with t-channel 
LF-boson. Numerically, this contribution amounts rel¬ 
atively to around — 2 % of Ap^ in the high invariant 
mass bin, and is insignificant elsewhere. Qualitatively, 
the charm and beauty asymmetry predictions are simi¬ 
lar. However, unlike the beauty predictions the mixed 
QCD-EW corrections to the asymmetry are positive, be¬ 
cause up-type and down-type quarks have opposite weak 
charges. This leads to a slightly increased asymmetry in 
both the low and high invariant mass bins. Another no¬ 
table difference is that the value of the asymmetry in the 
resonant bin is reduced for the charm-quark final state as 
the magnitude of the vector coupling of the Z boson to 
up-type quarks is approximately two times smaller than 
that to down-type quarks. 


To better quantify the differences between Ap^ 
and A|c, we present predictions for the ratio of these 
two asymmetries in Table EYl The differences in these 
asymmetry ratios due to corrections involving electro¬ 
magnetic/weak couplings range between 8 % to 30%, sug¬ 
gesting that an improvement in the experimental sys- 
tematics will be required to observe such effects — the 
systematic uncertainty in the central mass bin of the 
i/s = 7 TeV asymmetry measurement is about 30%. 


Another interesting feature of the bb and cc predictions 
which we wish to highlight is the relative size of symmet¬ 
ric cross section predictions, which is almost entirely due 
to QCD. Although the scale uncertainties are still sig¬ 
nificant at 0{a^), the ratio of charm and beauty cross 
sections is, for large values of the invariant mass of the 
heavy-quark pair, expected to be robust with respect to 
higher-order QCD corrections. We display our NLO and 
LO predictions for the ratio of symmetric cross 

sections in Table llVl ^ Given that this observable is the¬ 
oretically under control, and directly sensitive to charm- 


® At LO the scale uncertainties almost perfectly cancel between the 
numerator and denominator of so that the uncertainties 

shown in Table nv] for this case are not a reliable estimate of 
higher-order corrections. 
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tagging and beauty-tagging efficiencies, it should prove 
useful for validating these efficiencies experimentally for 
large values of rngg. 


CONCLUSIONS 

In this letter we have presented the results of an im¬ 
proved fixed-order computation for the forward-central 
asymmetry beauty-quark pair production, in¬ 

cluding a realistic evaluation of scale uncertainties. Our 
NLO predictions agree within uncertainties with the re¬ 
cent LHCb measurement of this quantity that has been 
performed with lfb“^ of a/s = 7TeV data. Anticipating 
improved LHCb measurements of the bb asymmetry, and 
a first measurement of its counterpart involving charm- 
quark pairs, we have additionally provided predictions 
for and Ap/, at ^/s = 13 TeV. 

Another interesting observable, which we have identi¬ 
fied, is the ratio of the differential cc and bb production 
cross sections at high invariant mass. Unlike many heavy 
quark-pair observables, this ratio is expected to be robust 
with respect to higher-order QCD corrections. This fea¬ 
ture can be used to validate charm-jet tagging efficiencies 
and mis-tag rates in future LHCb analyses. A very good 
understanding of these experimental issues is a prereq¬ 
uisite for any attempt to search for processes such as 
pp ^ h (—>■ cc) W, Z in the forward region. 
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